ECTP-2011-01 



(N 



Or 



The Hubble parameter in the early universe with viscous QCD matter and finite 

cosmological constant 

A. TawfikQ 

Egyptian Center for Theoretical Physics (ECTP), MTI University, Cairo-Egypt 

The evolution of a flat, isotropic and homogeneous universe is studied. The background geometry 
in the early phases of the universe is conjectured to be filled with causal bulk viscous fluid and 
dark energy. The energy density relations obtained from the assumption of covariant conservation 
of energy-momentum tensor of the background matter in the early universe are used to derive the 
basic equation for the Hubble parameter H. The viscous properties described by ultra-relativistic 
equations of state and bulk viscosity taken from recent heavy-ion collisions and lattice QCD calcu- 
. lations have been utilized to give an approximate solution of the field equations. The cosmological 

constant is conjectured to be related to the energy density of the vacuum. In this treatment, there 
is a clear evidence for singularity at vanishing cosmic time t indicating the dominant contribution 
, from the dark energy. The time evolution of H seems to last for much longer time than the ideal 

■ case, where both cosmological constant and viscosity coefficient are entirely vanishing. 
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vL, ■ I. INTRODUCTION 

Dark energy and dissipative processes are supposed to play a very important role in the evolution of the early 
CN ! universe and therefore would have essential astrophysical and cosmological consequences [l| . Much progress has been 
achieved in relativistic thermodynamics of dissipative fluids [H, Q. The widely used theoretical framework is the 
[ Israel and Stewart (IS) theory 4 L 5[, in which the causality is apparently conserved and the theory itself seems to 
. be thermodynamically stable 16 , 7 L The dynamics and evolution of the early universe can be described by causal 
bulk viscous thermodynamics [8l4ll|. But due to their complicated characteristics, very few exact solutions of the 
gravitational field equations are derived, especially in framework of full causal IS theory [I^ - [l6j . Also, it has been 
proposed that causal bulk viscous thermodynamics can model on a phenomcnological level matter creation in the early 
universe [HI, [l3j] . Recently, the effects of bulk viscosity on the early universe with vanishing cosmological constant A 
have been investigated 

Basically, the mathematical difficulties of the non-linear and non-homogeneous differential equations result in limiting 
most investigations of dissipative causal cosmologies to homogeneity and isotropy Fricdmann-Lcmaitre-Robertson- 
Walkcr (FLRW) symmetry The Einstein field equations for such models can be decoupled and therefore are 

■ reduced to an autonomous system of first order ordinary differential equations, which can be analyzed approximately 

To have a reliable description, the background matter in FLRW model is assumed to be filled with dark energy 
and an ultra-relativistic viscous QCD matter. The bulk viscosity of the QCD matter and its equation of state have 
been deduced from recent heavy-ion collisions experiments and the lattice QCD simulations. In the present work, we 
study such a sophisticated system of non-linear and non-homogeneous differential equations at finite bulk viscosity 
and cosmological constant. 

Recent RHIC results give a strong indication that hot dense matter is conjectured to be formed in the heavy-ion 
collisions experiments [181 ] . Such an experimental evidence seems to agree with the "new state of matter" as predicted in 
the lattice QCD simulations [l9j . However, the experimentally observed elliptic flow in peripheral heavy- ion collisions 
would give an indication that a thermalized collective QCD matter has been produced. In a addition to that, the 
success of ideal fluid dynamics in explaining several experimental data, e.g. transverse momentum spectra of identified 
particles, the elliptic flow poj shows that a nearly perfect fluid is likely to created and the quarks and gluons are likely 
to go through a relatively rapid equilibrium characterized with a thermalization time less than 1 fm/c pll ]. According 
to recent lattice QCD simulations [23 ], the bulk viscosity £ is not negligible near the QCD critical temperature T c . It 
has been shown that the bulk and shear viscosity at high temperature T and weak coupling ct s arc respectively given 
as £ ~ a 2 s T z /\aa~ 1 [23J. Such a behavior obviously reflects the fact that near T c QCD is far from being conformal. 
But at high T, QCD approaches conformal invariance, which can be indicated by low trace anomaly (p — 3p)/T 4 [24] ], 
where p and p are the energy and pressure density, respectively. 
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The present work is devoted to investigate the effects of dark energy and bulk viscosity on the evolution of the Hubble 
parameter H in the early universe, especially in the quark-gluon plasma (QGP) era, which is positioned within the 
range 0.2 < T < 10 GcV or 18.35 < t < 0.0073 GeV -1 . This limitation is likely compatible with the recent lattice 
QCD simulations and the experimental evidences that QGP seems to remain strongly correlated, i.e. non- viscous, up 
to several T c . Also, this sets on the upper and lower region of the validity of this treatment. We consider a background 
corresponding to a FLRW model filled with ultra-relativistic viscous QGP matter, whose bulk viscosity and equation 
of state have been deduced from recent heavy-ion collisions experiments and lattice QCD simulations. 

The present paper is organized as follows. The cosmological constant and the equation of state of viscous QGP are 
reviewed in section [IT] Section lllll is devoted to the evolution equations in viscous background QGP matter at finite 
cosmological constant. In section IPVl we present analytical solutions for the evolution equations. The results and 
conclusions are given in section [V] 



II. THE COSMOLOGICAL CONSTANT AND VISCOUS QGP EQUATION OF STATE 



A. The cosmological constant 

The cosmological constant A was a subject of various modifications. Originally, it has been introduced to fit 
astrophysical data [25f. The essential role of A is to allow for a static solution to Einstein's equations, Eq. (JSJ) , in a 
cosmological background filled with matter has been fundamentally modified, especially at the time of the discovery 
of the universe expansion (26j . 

The cosmic microwave background (CMB) observations indicate that the universe seems to have a negligible space 
curvature. Hereafter, we shall assume that the universe is spatially flat. Other astrophysical and cosmological data 
gathered over the last couple years apparently support the models with finite dark energy An example can be taken 
from the redshift of supernovae [27l ]. This apparently shows that the universe is currently undergoing an accelerated 
expansion. Furthermore, it seems that the flat universe contains both matter and cosmological constant A so that 
or for a stringent upper limit < 0.51 (at the 95% confidence level). On the other hand, there are many 
observations indicating that the non-relativistic matter contributes with ~ 1/3 of the critical density [28j . In other 
words, such observations would also mean that the universe contains ~ 2/3 of its total energy density in form of dark 
energy. The latter is a substance with negative effective pressure responsible for the current accelerated expansion 

Although, the original motivation for the cosmological constant A is not entirely clear [3l[ , it may be connected with 
scalar fields [32j (dilaton and moduli). Also, the interactions with the background strings may produce chaos [H| and 
hence A would gain small stochastic contributions (34[. A can be taken as a legitimate addition to the gravitational 
field equations, Eq. ([5|) for instance, or as a parameter to be constrained by observation, as mentioned before. The 
nature of dark energy itself is a subject of speculation. Apparently, it is very homogeneous but not very dense. It 
interacts through the gravitational force, exclusively. In order to determine the origin of A, all other terms in Eqs. ^ 
and (|10[) have to be measured to a sufficient precision. Such a way, one finds out that the A/3-tcrm has to be finite. 

From theoretical point-of-view, A can be interpreted as a measure of the energy density p of the vacuum which 
is supported by the existence of "zero point" energy in quantum mechanics based on particle-antiparticle pairs and 
confirmed through Casimir effect. The total energy density stems from various mostly unrelated sources. Each of 
them seems to come up with a contribution which is larger than the upper limits of today's A- value [35L |36| . In general 
relativity, all forms of p are eligible to gravitate through modifying the space. A is referring to p vac = A/8tt and a 
negative gravitation and therefore seems to explain the expansion nature of the non-empty universe. Furthermore, 
there is a recent suggestion that the velocity of light c and the exp ansion of the universe are two aspects of one single 



concept connecting space and time in the expanding universe 37 1. 

On the other hand, the role of cosmological constant in QGP era is conjectured to be played by the bag constant B 
(38j . It would be interesting to see what kind of long range effects (like interactions or correlations) B may produce. 
For the high density QGP phase, it seems to be eligible to add a cosmological term with a large cosmological constant 
of the order of the bag constant B to the covariantly conserved energy-momentum tensor of a viscous fluid. At low 
temperatures (< 10 MeV), a tiny cosmological constant is to be added to T Ml/ of an ideal mixture. An understanding 
of the transition between these two descriptions is, of course, completely in the dark. This is actually an aspect of the 
profound A-problem. On the other hand, such a description would require to have equations of state describing the 
two phases (hadrons and QGP), simultaneously. 

Based on various experimental evidences, it is assumed that B at very high temperatures likely disappears from the 
equation of state. At such high temperatures, the cosmological constant A has to be introduced, explicitly. Also, it 
should be noticed that the structural latent heat of QGP would lead to small value of A. As will be given in Eqs. ^ 
and pop . A is the only positive term in the field equations. Another stunning fact is that A remains constant, especially 
when all kinds of interactions between matter and radiation entirely disappear [39} . Its role apparently depends on 
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reconciling the age problem of the universe and the overall matter content [40J ■ The latter can be dominant or minor 
relative to A and therefore the evolution of different cosmological parameters would reflect the role of A. The presence 
of A apparently influences the expansion of the universe, i.e. its age. 

The reasons why A is to be included in QGP era are the long-range correlations and the constant self-interactions 
in QGP matter 1381 . besides the cosmological motivations. Recent results from RHIC approve the strong correlations 
in QGP phase |18j |. The equation of state of dark energy can simply be characterized by negative u, Eq. ([2]). In the 
case that oj = — 1, then 

* = (1) 

Recent observations of supernova, galaxy clusters, and CMB seem to confirm that uj has values slightly deviate from 
— 1 [4l|. In section that follows, the equation of state of viscous QGP is deduced. 



B. The equation of state of viscous QGP 



The equation of state, the temperature and the bulk viscosity of QGP, can be determined approximately at high 
temperatures [42j from recent lattice QCD calculations j4^|, as 

V = up, 

T = $p r , (2) 
£ = ap x , 

with u) = (7 - 1), 7 ~ 1.183, r ~ 0.213, (3 ~ 0.718, a = {9j 2 - 24 7 + 16)/((7 - l)w ) and cj ^ 0.5 - 1.5 GeV. In order 
to close the system of the cosmological equations, we still need an expression for the relaxation time r 0, 

T = ^ = ap x ~ 1 . (3) 
P 

Expressions ^ and ^ are known as barotropic equations. Such a dependence has a practical advantage. As we give 
below, the field equations, Eq. ^ and (JTUJ), relate the cosmological parameters, like Hubble parameter H and scale 
factor a, to the energy density p. Again, the expression ([2]) arc standard in analyzing the viscous cosmological models, 
whereas the equation for r is a simple procedure to ensure that the speed of viscous pulses does not exceed the speed 
of light. 

Using this equation of state sets the validity of this treatment. It deals with viscous QGP matter. In present work, 
we mean with early universe a well-defined era, where the background matter is characterized by QGP. 



III. EVOLUTION EQUATIONS IN VISCOUS BACKGROUND QGP MATTER 

We assume that the geometry of the early universe is flat and filled with a bulk viscous cosmological fluid, which 
can be described by a spatially flat FLRW type metric with c = 1 given by 

ds 2 = dt 2 - a 2 (t) [dr 2 + r 2 (dd 2 + sin 2 9d0 2 )] . (4) 

The Einstein gravitational field equations in Robertson- Walker metric with a positive cosmological constant read 



1 . 8ttG 



R ij - ^9ij R + 9i]h~ = Tij, (5) 



where A is seen as a free parameter. The energy-momentum tensor of bulk viscous cosmological fluid filling the 
background geometry is given by 

T( = {p+ Peff ) UiU 3 ~ p e ffSf . (6) 

Following the Einstein's prescription that the cosmological constant can be treated as an independent parameter so 
that the energy-momentum tensor of background vacuum is given as 

« ~ 8tt ' {> 
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where i and j take 0, 1, 2, 3 and p is the mass density. When taking Eq. ([7]) into account, Eq. ((5|) can be modified as 
follows. 

Ri S - ± 9ij R = 8ttG (T 4J - 9ij p vac ) , (8) 

where p vac = A/87T is the vacuum energy density. 

Generally, the bulk viscous effects can be described by means of an effective pressure p e tf which includes the 
thermodynamic pressure p and the bulk viscous pressure II Q. The four velocity it, is conjectured to satisfy the 
condition UiU 1 = 1. The particle and entropy fluxes characterize the background fluid, nearly entirely. Both are 
defined according to N l = nu % and S l = sN l — (tII 2 /2^T) u 1 , where n is the number density, s the specific entropy, 
T > the temperature, £ the bulk viscosity coefficient, and r > the relaxation coefficient for transient bulk viscous 
effect sometimes referred to as the relaxation time, Eq. ([3]). In AT = 4 SYM theory, the relaxation time r can be given 
by (2 — ln2)/(27rT) [HJ|. This has a very much little physical motivation. 

Taking into account a finite cosmological constant A, then in the comoving frame the energy momentum tensor has 
the components T® = p,T\ =T% = Tf = —p e ff- For the line element given by Eq. the Einstein field equations 
read 

3H 2 = 8nGp + A, (9) 
2,[H + H 2 ^ = -4nG (p + 3p eff ) + A, (10) 

where an over-dot denotes the derivative with respect to the cosmic time t and G is the gravitational constant. The 
Hubble parameter H is related to the scale factor a via H = a/a. In rest of this article, we take into consideration 
natural units, i.e. c = H = G = 1, for instance. The latter has been conventionally modified as, (47r/3)G = 1. 
Therefore, in order to obtain the expressions given in Ref. [ll|, A and A should be assigned to and 1, respectively. 

Assuming that the total matter content of the universe is conserved, TJ? . = 0, then the evolution of the energy 
density of cosmic matter seems to fulfill the following conservation law: 

p=-3H(p+ Peff ), (11) 

which simply can be obtained from Eqs. © and (TTU|) . Apparently, this relation entirely eliminates the cosmological 
constant A. 



A. Eckart relativistic viscous fluid 

The first attempts at creating a theory of relativistic fluids were those of Eckart and Landau and Lifshitz Q. 
These theories are now known to be pathological in several respects. Regardless of the choice of the equation of state, 
all equilibrium states in these theories are unstable and in addition signals may be propagated through the fluid at 
velocities exceeding the speed of light c. These problems arise due to the nature of the first order of this theory, that 
it considers only the first order deviations from the equilibrium leading to parabolic differential equations, hence to 
infinite speeds of propagation for heat flow and viscosity, in contradiction with the principle of causality. Conventional 
theory is thus applicable only to phenomena which arc quasi-stationary, i.e. slowly varying on space- and time-scales 
characterized by mean free path and mean collision time. 

The Eckart theory can be applied on modeling the cosmic background fluid as a continuum with a well-defined 
average 4- velocity field u a where u a u a = — 1. The vector number density n a = nu a can be estimated, when unbalanced 
creation/annihilation processes take place; n a a = 0. This apparently means that 

h + 3Hn = 0, (12) 

where the Hubble parameter H = Vu. In case of viscous fluid, the entropy current, 

S a = snu a , (13) 

is no longer conserved. The covariant form of second law of thermodynamics is > and the divergence of entropy 
current is given by TS" a = — 3HU. 

The evolution of the cosmological fluid is subject to the dynamical laws of particle number conservation Nl = 
and Gibbs' equation Tdp = d(p/n) + pd(l/n). In the following, we suppose that the energy-momentum tensor of 
the cosmological fluid is conserved, T\. x = 0. Then, based on Gibbs' equation, the covariant entropy current can be 
achieved by 



n = -3iH. 



(14) 
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This is a linear first-order relationship between the thermodynamical flux II and the corresponding force H. Substi- 
tuting with the system of equations given in Eq. ^ on the equation §§§ and (|10[) results in 

H = Una [ 8n ) H-- 1 H 2 + - 7 A. (15) 
Analytical solutions for this differential equations will be introduced in section HV Al 



B. Israel- Stewart relativistic viscous fluid 

A relativistic second-order theory has been developed by Israel and Stewart 0, Q , Hiscock and Lindblom Q through 
"extended" irreversible thermodynamics. In this model, the deviations from equilibrium (bulk stress, heat flow and 
shear stress) are treated as independent dynamical variables, resulting in 14 dynamical fluid variables to be determined. 
The causal thermodynamics and its role in general relativity are reviewed in Ref. Q- A general algebraic form for S a 
including a second-order term in the dissipative thermodynamical flux IT 0, Q reads 

S a = snu a + ftU 2 ^, (16) 

where ft is a proportional constant. 

For the evolution of the bulk viscous pressure, we adopt the causal evolution equation Q obtained in the simplest 
way (linear in II) to satisfy the iJ-theorem (i.e., for the entropy production to be non- negative, = IT 2 /£T > 
0, Hi ) • According to the causal relativistic IS theory, the evolution equation of the bulk viscous pressure reads Q 

rn + n=-3tff-irn^3F+I-i-^J. (17) 

In order to have a closed system from equations (J5J and (fTTj) . we have to take into consideration equations of state for 
p and T. 

With the use of Eqs. (fTTj) . © and ©, respectively, we obtain the following equation describing the cosmological 
evolution of the Hubble function H 

H I 3 B HH 2 I ( ^~ B l) | (£ + 6^7)A 2VCA 1 18tt 7 2 \ . _ 

2tt A y 2 3X A a{~A) x H A J 16tt 

277 VC , flB _27( 7 -2) g6+ /9 7 _ \ 3A VC + 37^ + 4,(3, - 8) ^ + 



2 4 ~ 3X a {-A) X A AA \A 'J 2 2 - 3A 7ra (-A) X A 8nA 

A 2 VC u2 , 4(B + 15,) + 3 7 (i? + 67r) A 



8 7 J 2 1 - 3 *7ra {-A) X A H ~ + 24irA A H 



7 1\ A 2 VC 



6 J 2 1 - 3A ttq (-A)M 

where ^4 = — 77 2 + A, B = 1 + r and C = 7r 2A . The analytical solutions are given in section HVBl 



0, (18) 



IV. ANALYTICAL SOLUTIONS AT FINITE COSMOLOGICAL CONSTANT 
A. Analytical solutions in Eckart relativistic viscous fluid 



At arbitrary values of the parameters A and A, the expression (| 1 5[) apparently turns to be non-integrable. Let us 
first check the default values A = and A = 1. The solutions simply read [ll[ 



3 1 



7 + 3aH In ( 3a — — 



(19) 



On the other hand, when assuming that A = 1 + e, where e is a small positive value, then Eq. (|15p can be reduced to 

(20) 



-aH 3 --jH 2 --aAH+-^A, 
2 2' 2 2 ' ' 
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which can be solved as 



1 



3 7 2 - 9a 2 A 



2 , y\j — arctanh 



(21) 



The graphical representations of Eqs. (fl9|) and (f2~Tj) are given in Fig. [3] The first one seems to reflect the treatment 
of Rcf. [ljj, where the dark energy, A, has been excluded. 

Although the mathematical difficulties arising with utilizing IS theory in the early universe, its application is simply 
unavoidable. It is derived by the real concern about the causality constrains and non-consistency with the second law 
of thermodynamics. 



B. Analytical solutions in Israel-Stewart relativistic viscous fluid 



We introduce the transformation u = H, so that Eq. (|18p is transformed into a first order ordinary differential 
equation, 



du 3 B 

u 1 Hu 

dH 2tt A 



, 3(8tt - B-y) (B + 6tt 7 )A 2s/CA 1 18tt 7 , 

H 



2 3A 



277 VC tt6 27( 7 - 2)^ 5 



24-sAq, (_a) x A 



H 



4.4 



.4 
3A 



a(-A) x H A 



2 3A 
16^ 



-Hu 
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4 7 ~ J 2 2 - 3A 7ra (—A) X A 



C g 4 + 37^ + 4^(37-8)^3 
8ttA 



C 



2 1 - 3A ttq [~A) X A 



2A-kA 
7 1\ A 2 



C 



6 y 2 1 - 3A ttq (— A) A A 



Using the new variable 



we can rewrite Eq. ([22)) in the form 



fl = uE=uexp(-J ^-~J HdH ) 



dtt 
dH 



= 0. (22) 



(23) 



where 



Fi(H) 
F (H) 



2- i+3X (-A)£ 

7T 

(-A)- 1+ £- x 
48na 



3 8 



in - B 7 ) 



(5 + 6tt7)A 2,4^ 18H 2 irj 



H. 



A a{-A) x H A 

X VC (-81iJ 6 7r 7 + 9H 4 (-4 + 9 7 )A + (-4 + 3H 2 (8 - 9 7 ) + 3 7 ) A 2 



(24) 



(25) 



-2{-A) x Ha (162H 4 tt(-2 + 7) - 9H 2 (Bj + 4tt(-8 + 3 7 ))A - (4B + 60tt + 3{B + 67r) 7 )A 2 )] . (26) 
By introducing a new independent variable z = J F\ (H) dH, we obtain 



where g(z) — F /Fi can be approximated as a simple function depending on the new variable z as follows. 



(27) 



9(H) 



z(H) 



-A) 



_b_ 8 q VC (-81if 6 7T7 + 9i? 4 (9 7 - 4)A + (3iJ 2 (8 - 9 7 ) + 3 7 - 4) A 2 ) 
-3 2 1 + 3 iA 2 VC + 3(-A)iHa (3A(8ir - B7) + 8<? (BA + 6.7,4)) 
_b 2(-A)iHa (I62i? 4 7r(7 - 2) - 9H 2 (Bj + 4.(37 - 8))A - (45 + 60tt + 3(5 + 67r) 7 )A 2 ) 



->3A-5 



(-A) = 



.4 



-3 2 1 + 3 iA 2 VC + 3(- A)iHa (3A(8. - B7) + 8<- 
-9. 7 # 3 + ( J^Stt - Bj) + (B + 6tt 7 )a) ff - ^ ( _ A)A _ 2 



(.BA + 6.7A)) 
^ 5. 



(28) 



(29) 
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Fig. 1: In left panel, g(H) = Fq/F\ is given in dependence on H, whereas in right panel z{H) is depicted. 




Apparently, replacing H with z is a non-trivial task. To this goal, let us first depict the real parts of both functions 
g(H) and z(H), Fig. [T] Obviously, the two functions have various ingredients. Now, we may try the parametric 
plotting, as in Fig. O It is obvious that the dependence of g on z can approximately be characterized as linear. 



where 



g(z) « (f>z, 

-0.07 ± 0.02. Then, from the definition of SI, we simply have 



n = (3H 2 -A) 



-B/4tv 



H. 



(30) 



(31) 



Now, we reduce the last expression to the canonical equation of the Abel type. In doing this, we use the relation 
SI = zjV . Then from Eqs. (|3"Tj) and (|29[) . we obtain a first order differential equation for H. 



->3A-5 



VH 



H 



-9Tr-fH 3 + ( ^A(8tt - Bj) + A(B + 6tt 7 ) j H - 



4VC 



2-A 



This equation has the solution 
V 



t 



'.KM 



3a(B-f - 8tt) + 8 A ( 12VC(X - 2) \n(H) H — a(B + 6nj) 



32tt 
If 



348tt 



(32) 



(33) 



where Af = (247T — 'SBj + 8 X (B + 6ttj)) 2 and V is taken as a free parameter. We can assign any real value to V . For 
the results presented in this work, we used a negative value. This negative sign is necessarily to overcome the sign 
from the integral limits. 



V. RESULTS AND CONCLUSIONS 



In the present work, we have considered the evolution of a full causal bulk viscous flat, isotropic and homogeneous 
universe with finite bulk viscosity parameter (£ ^ 0) and and cosmological constant (A ^ 0). The equation of state of 
QGP is taken from recent lattice QCD simulations and heavy-ion collisions. 

In Fig. [31 the cosmic time t and Hubble parameter H in cosmic background filled with the Eckart fluid arc depicted. 
In this treatment, t is given in GeV -1 whereas H in GcV. The solid curve represents the analytical solution, Eq. (12T1) . 
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0.35 
0.30 
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0.20 
0.15 
0.10 
0.05 




H 



8 



10 



Fig. 3: Cosmic time t vs. Hubble parameter H in cosmic background filled with Eckart fluid. The solid curve represents the 
solution in Eq. (|2ip. The solution at vanishing viscous and A is given by the dotted curve. The dashed curve represents Eq. 
(fT9)l . where A = and A = 1, i.e. at vanishing dark energy. 



Using equation of state of an ideal and non- viscous gas and a vanishing cosmological constant A result in the dotted 
curve. The dashed curve gives that results when A is vanishing and q = 1. This has been introduced in Ref. (Tlj . 
In all curves, the singularities are present. On the other hand, the time evolution shows essential differences. Eckart 
fluid enables H to decay much faster (solid and dashed curves), especially at small t. This behavior is flipped at large 
t-values. This might reflect the effects of the approximation in A = 1 + e, where e is taken to be very small. 

In Fig. 21 the analytical solution given in Eq. (|33[) is graphically illustrated. To have a comparison with other 
cases, we plot both non-viscous and viscous solutions. The dashed curve gives the latter case, where the cosmological 
constant is vanishing fill ]. 



t 



In 



1 



ar H 



1 



1 



+ 7 



(34) 



It is obvious that the causal bulk viscous universe described by this solution starts its evolution from an initial non- 
singular state with a non-zero initial value of H(t). The ideal case is characterized by vanishing viscosity coefficient £ 
and cosmological constant A, 



t 



3-yH 



(35) 



The singularity is apparently present. 

So far we conclude that the bulk viscosity seems to play an important role in the evolution of the early universe. 
Despite of the simplicity of our model, it shows that a better understanding of the dynamics of our universe is 
only accessible, if we use reliable equation of state in order to characterize the matter filling the cosmic background 
geometry. Also, the cosmological constant seems to be essential. The motivation of finite A is based on theoretical 
needs and astrophysical and cosmological observations [25M30l . It is very essential to describe the early stages of 
universe, Eq. ((TJ) . 

Assuming that A vanishes (Fig. 2] and Eq. (fM]) ). the solution shows an absence of singularity near t = 0. This 
was the main result of Ref. fill ] . Whether it has a physical meaning, we like to refer to the restriction to the QCD 
era of early universe. It should not be understood that it is applicable either in former or later eras. By way of 
precaution, the QCD era can be limited to the temperatures in the range 0.2 < T < 10 GcV or time in the range 
18.35 < t < 0.0073 GeV" 1 . Apparently, this limitation is likely compatible with the recent lattice QCD simulations 
and the experimental evidences that QGP seems to remain strongly correlated, i.e. non- viscous, up to several T c . 



9 




Fig. 4: Cosmic time t vs. Hubble parameter H in cosmic background filled with IS fluid. The solid curve represents the solution 
in Eq. (|33[) . The dashed curve depicts the results of Ref. [lj, where £ 7^ but A = 0. The solution at vanishing viscous and A 
is given by the dotted curve. Obviously, the negative values are non-physical. 



In the present work, it is assumed that the universe is spatially flat and the background geometry is filled with dark 
energy and QCD matter (QGP) with a finite viscosity coefficient. The results seem to support the singularity near 
t — > 0. The resulting universe is obviously characterized by almost the same behavior of the ideal case. There is two 
essential differences. First, the solution given in Eq. (|33jl depends on free parameter V . The solid curve can be moved 
up and down depending on V . Its shape is not depending on it. Second, the universe in which dark energy is taken 
into account seems to live much longer that the ideal case. 

The validity of our treatment depends on the validity of the equations of states, Eq. [21 which we have deduced from 
the lattice QCD simulations at temperatures larg er than T c w 0.19 GeV. Below T c , as the universe cooled down, not 
only the degrees of freedom suddenly increase [45| but also the equations of state turn out to be the ones characterizing 
the hadronic matter. Such a phase transition - from QGP to hadronic matter - would characterize one end of the 
validity of our treatment. The other limitation is the very high temperature (energy), at which the strong coupling 
a s entirely vanishes. 
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